Anomalous Hall effect in a two-dimensional electron gas with spin-orbit interaction 
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We discuss the mechanism of anomalous Hall effect related to the contribution of electron states 
below the Fermi surface (induced by the Berry phase in momentum space). Our main calculations 
are made within a model of two-dimensional electron gas with spin-orbit interaction of the Rashba 
type, taking into account the scattering from impurities. We demonstrate that such an "intrinsic" 
mechanism can dominate but there is a competition with the impurity-scattering mechanism, related 
to the contribution of states in the vicinity of Fermi surface. We also show that the contribution to 
the Hall conductivity from electron states close to the Fermi surface has the intrinsic properties as 
well. 

PACS numbers: 73.20.Fz; 72.15.Rn; 72.10.Fk 



I. INTRODUCTION 

The measurement of anomalous Hall effect 1 (AHE) in 
ferromagnetic metals and semiconductors is a powerful 
tool to characterize the magnetization of thin films and 
mesoscopic objects. That is why it attracted recently a 
lot of attention in the spintronics community.^ On the 
other hand, the theoretical aspects and possible physical 
mechanisms of AHE are still not quite well understood, 
and the great efforts are directed now toward the physical 
explanation of AHE in magnetic semiconductors, pyro- 
clore crystals, and magnetic compounds with inhomoge- 
neous magnetization. The state of art in the theory of 
AHE has been recently reviewed by Sinova et al& 

Relatively well studied interpretations of AHE include 
the skew scattering^ and side-jump^ mechanisms, which 
have been recently revised and physically clarified by 
Crepieux et al£& These mechanisms are related to the 
scattering from impurities, with the spin-orbit (SO) in- 
teraction playing the leading role, so that the spin-flip 
interaction is included into the amplitude of scattering 
from impurities. 

It was shown recently that the nonzero spin chiral- 
ity in some inhomogeneous ferromagnets also leads to 
the nonvanishing AHE. 910,11 i*i2ii 3 . Such a chirality mech- 
anism of AHE is related to the Berry phased acquired 
by the wave function of electrons moving adiabatically in 
the inhomogeneous magnetization field. In principle, the 
SO interaction is not a necessary element of this mech- 
anism. However, in earlier presented theoretical models, 
the nonvanishing average AHE was obtained only in the 
presence of SO interaction. Otherwise, the spatially av- 
eraged off-diagonal conductivity was found to be zero. 
A model of magnetic nanostructure without SO inter- 
action, leading to the Berry phase induced AHEji£ has 
been called the topological mechanism of AHE. 

Now the main interest is directed toward another "in- 
trinsic" mechanism of AHE recently proposed by Mac- 
Donald et a?»iS*ii and Nagaosa et alm& This mechanism 
uses the SO interaction included into the Hamiltonian of 



electrons in the crystal lattice. The key element of this 
theory is the Berry phase in the momentum space, so 
that a nontrivial topology of the electron energy bands is 
necessary for the intrinsic mechanism. A similar idea of 
the calculations of Hall conductivity has been used long 
ago in the context of quantum Hall effect piiLSfli 

The Berry-phase-induced intrinsic mechanism of AHE 
does not take into account scattering from impurities. It 
was suggested before that the impurities can be totally 
neglected in case of intrinsic SO interaction because the 
effect of impurities should be necessarily extrinsic. Corre- 
spondingly, the calculations of AHE have been performed 
in a zero density-of-impurity limit, Ni — > 0^£*i& which 
was taken before the static limit of lu — ► in the Kubo 
formula for conductivity o~ij(uj). 

In this paper we present the result of calculations tak- 
ing into account the effect of impurity scattering. We 
demonstrate that there exists another intrinsic contribu- 
tion to AHE related to the electron states near the Fermi 
surface. This contribution is not geometric, and we show 
that it can be of the same order of magnitude as the 
Berry-phase-induced one, which has been calculated be- 
fore. To calculate properly the non-geometric contribu- 
tion, we should take into account scattering from impuri- 
ties leading to the finite relaxation time of electron states 
at the Fermi surface. Accordingly, we take the limit of 
u) — > before Ni — » 0, which ensures the transition be- 
tween disordered and clean to be continuous. At the same 
time, the non-geometric contribution is intrinsic, and it 
does not vanish if we finally take the limit of Ni — > 0. 

In our calculations, we use a model of two-dimensional 
electron gas (2DEG) with SO interaction in the form of 
Rashba termini We show that, in accordance with the 
well-known result of Stfeda^ii there are two terms in 
the off-diagonal conductivity, one of which, a 1 , is due to 
the electron states near the Fermi energy (it corresponds 
to the non-geometric contribution), and the other one, 
cr^y , is related to the contribution of all occupied electron 
states below the Fermi energy (identified as the Berry- 
phase-induced intrinsic mechanism, which has been cal- 
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culated before). 

In agreement with other works on different models of 
the electron energy spectrum, our calculation confirm 
that for the 2DEG with Rashba SO interaction, there is 
a non-vanishing contribution <r~.y from the states below 
the Fermi surface. This contribution is absent in the case 
of the side-jump and skew scattering of electrons within 
the model of a parabolic energy band. However, we em- 
phasize that the contribution to a xy from the vicinity of 
the Fermi surface (cr^ ) can not be totally neglected. 

Recently, the role of impurity scattering in AHE has 
been analyzed within the double-exchange model with 
non-coplanar spin configuration^ It was shown that two 
different mechanisms can contribute to the AHE, one of 
which is related to the spin texture in the real space and 
the other to a momentum-space skyrmion-density at the 
Fermi level. In the latter case, the impurity scattering 
was relevant. This conclusion is in agreement with our 
results on a different model. 



II. MODEL AND THE KUBO FORMULA 

We consider a model of 2DEG in the x—y plane under a 
homogeneous magnetization field Mo directed along the 
axis z, with the SO interaction described by the Rashba 
term^Ai The Hamiltonian of this model reads 



H = £k + a (a x k y - a y k x ) - Ma z 



(1) 



where £k = k 2 /2m, a is the coupling constant of SO 
interaction, M = g^sM®, Mq is the magnitude of mag- 
netization, and we take units with h — 1. Due to the 
SO interaction, the components of velocity operators, 
Vi = dH/dki, are matrices in the spin space 



k 



V 

m 



(2) 



We also include into consideration the scattering from 
impurities, described by a disorder potential V(r). We 
assume that the disorder potential is short-range and 
weak, and we will treat it in the Born approximation 
of the impurity scattering. Thus, the Hamiltonian of our 
model is H = Ho + V(r), and any physical variables in- 
cluding the Hall conductivity should be calculated with a 
corresponding averaging over the random potential V(r). 

It should be noted that such a model does not take 
into account the SO interaction in the impurity poten- 
tial V(r). It is known that the SO interaction in the 
impurity scattering produces the AHE via well-studied 
side-jump and skew scattering mechanisms (see Ref. 
and references therein). Here we concentrate on other 
possible mechanisms induced by the SO interaction in 
the potential of the crystal lattice. 

The conductivity tensor can be calculated using the 
general Kubo formalism 



where G(e) = (e — H)~ is the operator Green function 
of electrons described by the Hamiltonian H, which in- 
cludes the disorder, di is the corresponding velocity op- 
erator, the trace goes over any eigenstates in the space 
of operator H (e.g., eigenstates of the Hamiltonian H 
itself), and (...) means the disorder average. 

It was first shown by Stfeda™ and later in context of 
the AHE by Crepieux et ali that a general formula for 
the static conductivity ov,(w = 0) can be presented in 
a form of two different terms, a xy — u xy + a~,y , one of 
which, a^y , stems from the contribution of electrons at 
the Fermi surface 



j e 2 fde ( df(e) 



G R (e) - G A (e) 



x Vj G A (s) - ^ G R (e) vj \G R {e) - G A (e)J \ , (4) 

and the second one, a xy , formally contains the contribu- 
tion of all filled states below the Fermi energy 



11 2 "J 2vr 
„ Af . A dG A (e) 



r>GUc)>\ — +v t G R (e)i 3 dGR{£) 



ds 



dG R (e 



r.G": i). (5) 



In this paper, for the model of Eq. (1) with the im- 
purity disorder, using a different way of calculation, we 
come to a similar result in accordance with the Stfeda 
formula. We calculate the off-diagonal conductivity a xy 
in the loop approximation, changing the disorder average 
in Eq. (3) as (v x Gv y G) — > T x (G) v y (G), where T x is the 
renormalized velocity related to the vertex corrections^ 
In this approximation we do not take into account the 
localization corrections f 2 * which are vanishingly small in 
the limit of small concentration of impurities. The ver- 
tex corrections are known to vanish in the case of short- 
range impurity potential and simple parabolic energy 
spectrum. As we show in Sec. 3, in the case of a complex 
band structure, the vertex corrections do not vanish even 
for the short-range impurity potential. 

To simplify the notations, in the following we omit the 
angular brackets for (G), and we denote the disorder- 
averaged Green function as Gk(e). Here we use the 
momentum representation for Gk(e) since the disorder- 
averaged Green function is diagonal in this representa- 
tion. 

Hence, we will calculate the off-diagonal conductivity 
from 



Ca;y(w) 



2tt (2tt) 2 



(3) 



Tr / ——^%{e,u)Gk{e + u)v y Gk{e), 

where the trace runs over the spin states. 

In the following calculation of the static off-diagonal 
conductivity a xy (uj = 0), we will be interested in the case 
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when the density of impurities Ni is small but finite. It 
also includes the "clean limit" of iVj — > 0, which can be 
physically realized in samples with vanishingly small con- 
centration of impurities and defects. Correspondingly, 
when we calculate the static conductivity tensor of a 
clean sample using Eq. (6), we take the limit of u> — > 
before the limit of Ni — > 0. 

It should be noted that in all previous considerations 
of the intrinsic mechanism of AHE, the impurities were 
totally neglected under an assumption that the effect of 
impurities is always " extrinsic" . Here we emphasize that 
this was essentially wrong. Performing the calculations, 
which include the scattering from impurities, here we find 
an additional contribution, which is also intrinsic, i.e., it 
does not vanish in the limit of Ni — + 0. This result is quite 
similar to the known side-jump mechanism of AHEi As 
we will see, the additional contribution to AHE comes 
from the states at the Fermi surface, for which the finite 
relaxation time due to impurity scattering is important. 
If we totally neglect the scattering from impurities, we 
miss this term. The real structures always have some 
impurities or defects. Therefore, we can justify taking 
the limit of oj — > before Ni — > as a physical realization 
of the clean system as a system with very small but still 
nonzero density of impurities. 

In the absence of impurity scattering, the electron 
Green function can be found using Hamiltonian (1) 



Gi(e) = 



£ - £k + M + a (k y a x - k x a y ) - Ma 
e — -Ek,+ + [i + iS sign e 
1 



e — £k,- + fi + iS sign e 



(7) 



where fi is the chemical potential and the electron energy 
spectrum consists of two branches, which we label as 
"+" and "-" corresponding mostly to the spin up and 
down electrons, respectively (even though they contain 
an admixture of the opposite spin due to the SO interac- 
tion), 



Ek,± =£kTA(fc), 



(8) 



1/2 



and we denote A(fc) = (M 2 + a 2 k 2 ) 

We consider a general case when the chemical poten- 
tial \x can be situated in both spin up and down sub- 
bands, corresponding to fi > M , as shown schematically 
in Fig. 1. When only the spin up subband is filled with 
electrons, —M < /x < M, only the contribution of the 
filled subband, -Ek,+ , should be hold. 

Due to the scattering from impurities and defects, the 
Green function (7) is modified. It is important to in- 
clude the effect of scattering for the correct evaluation of 
the contribution to the off-diagonal conductivity from the 
Fermi surface. For simplicity, we consider the model of 
disorder created by weak short-range scatterers, which 
can be treated in the Born approximation. The corre- 
sponding self energy of electrons is calculated as 



Si (e) = Ni 



Gi{e), 



(9) 




FIG. 1: Energy spectrum of electrons in a two-dimensional 
ferromagnet with Rashba SO interaction (schematically). 



where Vo is the Fourier transform of the impurity po- 
tential, and Ni is the impurity density. The result of 
calculation using Eqs. (7) and (9) can be presented as 



(e) = - 



i sign e 



1 1 

T T 



where 



1 



(10) 



(ii) 



(12) 



A± = X(kp t ±), and 



7)1 

2^ 



ma 



A + 



9{n - M) 



2tt 



1 - 



A_ 



(13) 



are the densities of states at the Fermi surfaces of two 
different subbands for /x > — M. Here 0{x) is the Heav- 
iside step function, and /cf,± are the Fermi momenta of 
the majority and minority electrons, respectively. 

Taking into account the self-energy correction (10), we 
find the Green function of the electron system with im- 
purities 



G k (e) 
where 

r = 



+ iT — £k + /x + a(k y tj x — k x (T y ) — ff z (M + iT) 



k.H 



i sign £ I 2t.| 
1 



£ — i?k,- + n + i signe/2r_ 



(14) 



sign £ 
2r 



signe 
2f 



1 



1 1 M 

T T X± 



(15) 



The values of t + and r_ determined by Eq. (15) are the 
relaxation times of electrons in different subbands. 
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III. VERTEX CORRECTION 



The equation for the renormalized vertex T(e, uj) can 
be presented using the diagrams with the impurity lad- 
der included into the vertex part^ For a short-range 
impurity potential, this equation has the following form 



Tifaw) =v i + N i Vg 



d 2 k 

(2^)2 



G k (e)7i(e,w)G k (e + a;). 

(16) 

In the limit of oj — > 0, the integral is not zero only at the 
Fermi surface, i.e., for energies e <C r^j. We assume that 
the density of impurities is low, which corresponds to the 
large relaxation times r-j- ^. Thus, we calculate the vertex 
renormalization for e = 0, otherwise we take T x — v x . 

Denoting T x = T x (e = 0, u> — » 0), we look for a solution 
of Eq. (16) in the form 



ca v 



(17) 



with some coefficients a, b, and c. After substituting this 
expression into (16) and using Eq. (14), we find that in 
the limit of low impurity density, a = 1/m and 6 = 0, 
and the constant c is defined by 



1 



NiV Q 2 



d 2 k (p-e k ) 2 



M 2 



1 + 2Ni Vr 



(2tt) 2 D+ D- 

d 2 k e k (fi - E k ) 



D+ D. 



(18) 



where D± = (ji- E K+ ± i/2r+) (/x - £ k - ± V 2r -)- 

Using Eqs. (2), (6), (17) and (18), we conclude that 
the vertex renormalization in Eq. (6) for a xy results in 
a substitution of the velocity v x by T x — k x /m — a*a y 
where 

J (2tt) 2 D+D_ 



1 - Ni Vt 



d 2 k (/!-£ k ) 2 -M 2 



(27T) 



D+ D- 



(19) 



We should emphasize that the vertex renormalization 
refers only to the states near the Fermi surface. 

The integrals in Eq. (19) can be calculated by trans- 
forming them to integrals over e k , i.e., J d 2 k (2n)~ 2 ... = 
J H)(ek) <fek—, where Vq{e) — (m/2n) 8(e). Finally, we 
find 



2vN t V 2 
(Si - £ 2 ) 2 



vo(£i) 



Si (/x - £{) 



2^0(^2) 



£2 (/x - £2) 



1 - 



nNi V 2 



(£1 



£2? 
-v$(£i 



vq{£ 



r 2 

[iv-£i?- 



M 2 



[(/x-g 2 ) 2 -M 2 ] 



(20) 
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FIG. 2: Renormalized constant of SO interaction a* in the 
vertex T x at the Fermi surface for different locations of the 
Fermi level. 



where 

51,2 = M + mo? i (to 2 a 4 + 2/ima 2 
and 



M 



(21) 



r lf2 = t 



H - £1,2 - (M 2 + 2ma 2 £ K 



.1/2 



2t+ (£ x - £ 2 ) 



F ■ 



fi - £12 + (M 2 + 2m.a 2 f i,. 



,1/2 



2t_ - £ 2 ) 



(22) 



If /i > M (Fermi level is located within two subbands 
like shown in Fig. 1) and in the limit of a — > 0, we can 
find from (22) that a* /a -> 0. 

In the case of — M < /x < M (only the lowest energy 
subband is partly filled with electrons) and a — > 0, we 
obtain 



a* _ ! /x + M 
~a~ ~ 2Af 



(23) 



which gives us a* /a — * 1 for /x — > — M, and a* /a — > 
for /i — > M. The dependence of a* /a on a for different 
values of /x is shown in Fig. 2. 
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IV. ANOMALOUS HALL EFFECT AND 
TOPOLOGY OF THE ENERGY BANDS 



The integration over energy e in Eq. (6) with Green 
function (14) leads to the separation of a xy into 
two parts, a xy = a xy + <r xy , corresponding to the 
above-mentioned contribution of energy states at the 
Fermi surface and the states below the Fermi energy, 
respectively^*^ 

First we calculate the contribution from the states be- 
low the Fermi energy, a^ y . Using (2), (6) and (14), we 
find a^y as an integral over momentum with the Fermi- 
Dirac functions f(E^±). Thus, it accounts for the con- 
tribution of all states with £k.+ < M an d E^ _ < ji: 



-Ma, 



ka- 
rri 



2ie 2 



a a. 



Tr 



d 2 h f(E K+ )-f(E K . 



(2tt) 2 (£ k ,+ - E K _y 



(E u 



ak y cr x 



+ aa x ) (—u + E K+ - e k 



ak x a y 



ak y o x 



In the static limit of lo —. 
trace, we obtain from (24) 



-ak x a y - Ma z ). (24) 
0, and after calculating the 



<r£l = -Ae 2 Ma 2 



d 2 k f(E K+ )-f(E K _) 



(2tt) 2 (E k+ - S k ,_)3 



(25) 



Using Eq. (8) we can calculate the integral over momen- 
tum and obtain finally 



xy Air 



1 



A 



(26) 



In the limit of weak SO interaction, akp.± <CM,we get 
from (26) 



jl 



9 1 

e ma 



2tt M 



6(M - n) 



M 



2M 



0{/j - M) 



(27) 



The expression for a x J y can be also presented in a differ- 
ent form demonstrating the topological character of this 
contributioniSi Let us introduce a 3D unit vector n(k) 
at each point of two-dimensional momentum plane 



n(k) = 



M 



( ak y ak x 
\X(k) ' ~Mk) ' ~ A(fc) 



(28) 



By using the n- field (28), we parameterize the manifold 
of 2 x 2 Hermitian matrices corresponding to Hamiltonian 
(1), since Ho = e k + X(k) <j ■ n(k). At k = 0, the vec- 
tor n is perpendicular to the k-plane, whereas for large 
|k| M/a, it lies in the k-plane and is oriented perpen- 
dicular to the momentum k, like schematically presented 
in Fig. 3. The dependence n(k) is a mapping of the k- 
plane to the unit sphere S2. As we can see from Fig. 3, 
the total k-plane maps onto the lower hemisphere of S^. 



/ / i 
i i \ 



T7 



-x. 

•V X V V 
V \ \ \ \ 



/ 



/ 
/ 



1 
/ 
/ 



FIG. 3: The field n(k) in the momentum space. 



Using (24) and (28) we find that in terms of the n-field, 
the contribution a x ! y can be written as 



'xy 



d 2 k 



[f(E K+ ) - /(£ k ,_)] 



where e Q /3 7 is the unit antisymmetric tensor. The integral 



2 J (2n) 2 



-f(E K -)} 



^^^^-dk-y (30) 

is the spherical angle on S2 enclosed by two contours L + 
and L-, where L± are the mappings of the Fermi surfaces 
(circles) -E kjr ,+ = (J- and -Efcjr,- = /x to the sphere S2, 
respectively. 

Notice that for — M < /i < M, and in the case of 
M/a < k F ,+ < 2ma(l + M/ma 2 ) 1 / 2 , we obtain fl ~ 2ir, 
corresponding to the mapping of the k-plane to the whole 
unit hemisphere. This is possible only in the case of 
large SO interaction or very small magnetization, when 
M < ma 2 . 

There is also a contribution to the Hall conductivity 
from the vicinity of Fermi surface. The calculation of 
corresponding term using Eqs. (6), (14) and (15) gives 
us 



a\. y = —e 2 Maa* 



d 2 k 



(27T) 2 (E K 



df(E K+ ) 

de 



1 - 



df(E K _) 
de 



l- T -± 



(31) 



Here the presence of factor (—df/de) restricts (31) to 
the integral over Fermi surface (we consider only the low- 
temperature limit of T/(/i + M) <C 1). Calculating the 
integral over momentum (31), we find 



e 2 aa*M 



■xy 



1-I±) 

rj A 2 (fc Fj+ ) 



+ 1- 



A 2 (fc Fi _) 



(32) 
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FIG. 4: ai y in units of e 2 /h as a function of SO coupling FIG. 5: al y in units of e 2 /h as a function of a for different /i. 
constant a for different filling of the energy bands. 



Using (11), (12) and (15), we can present a* depend- 
ing only on the parameters of the energy spectrum and 
density of states at the Fermi level 



e 2 aa" 



M 2 (£+-£_) 



A+[i/+ + i/_+M(£ 4 .-£_)] 



/_-M(£ + -£-)] 



,. y does not depend on impurities at all but is 



where £± = v±/\±. As we see from Eq. (33), the contri- 
bution of a], 

an intrinsic property of the crystal 

It follows from (33) that in the limit of large magneti- 
zation (weak SO interaction), M 3> akF,+ , the contribu- 
tion of this term is 



e maa 

47T 



M 



B{M-n)'- 



M 



40(jU-M) 



ma 



(34)" 

The dependence of two contributions to the Hall con- 
ductivity on the magnitude of coupling constant a for 
different positions of the Fermi level are presented in 
Figs. 4 and 5. We calculated them using the parameters 
M = 0.01 eV and m = m . Note that the contribution 
of a^y does not depend on the impurity density but only 



on the ratio between the relaxation times for majority 
and minority carriers. The figures demonstrate that a^y 
and a^y can be of the same order of magnitude, and for 
/i < M, the effect of a^. y can dominate. The total Hall 
conductivity is presented in Fig. 6. It shows that the 
AHE can change the sign with the density of carriers. 

^22^ The magnitude of SO coupling constant a can be 
roughly estimated if we take the SO splitting in a semi- 
conductor like GaAs, AEso — ak — 0.5 meV corre- 
sponding to the Fermi momentum k = (2-ttNs) 1 ^ 2 for 
the carrier density N s = 10 11 cm" 2 as a characteristic 
value for GaAs/GaAlAs heterostructures; 2 ^ It gives us 
a ~ 6.3 x 10~ 10 eV-cm. The experimentally obtained 
magnitude of a in InGaAs/InAlAs heterostructures is 
about 10 -9 eV-cmml 



V. SEPARATION OF THE CONTRIBUTIONS 
TO HALL CONDUCTIVITY: GENERAL CASE 



Let us consider now a general case of the energy spec- 
trum described by a Hamiltonian iJo, which is a matrix 
in a certain basis of wave functions. The off-diagonal 
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FIG. 6: Total Hall conductance a xy = <Jx V + crly in units of 
e 2 /h as a function of a for different fx. 



conductivity can be calculated using the Kubo formula 
a xy (u) = ^- J ^ ^ i v x)nm 

k nn'mm' 

xGk TOm '(e + w) (v y ) m i n > G k „'„(e), (35) 

where we can take the basis of eigenf unctions of the 
Hamiltonian H and assume a proper renormalization 
of velocities due to the impurity corrections. In this rep- 
resentation and in the absence of impurities, the Green 
function G£(e) is diagonal. However, the scattering from 
impurities induces nonzero off-diagonal elements, as we 
can see from the above-considered model (1) leading to 
the self energy with nonvanishing off-diagonal matrix el- 
ements after diagonalization of the Hamiltonian Hq. 

The off-diagonal matrix elements of the Green function 
can be formally presented in the following form resulting 
from an expansion over off-diagonal self-energy matrix 

Y](°//) 



where f„ 



_ zsigne P nn , (k) 

fnn' Dk(e) 



i n'n' (^) 

n ^ n' , 



(36) 



is the relaxation time associated with 
impurity-induced interband transitions, P nra '(k) is a co- 



efficient determined by the specific form of Hamiltonian, 
D k(s) = Iln ( £ - E kn +isigne/2r n ), E kn is the disper- 
sion of electrons in the n-th subband, and r„ is the cor- 
responding relaxation time, which includes all intra- and 
interband transitions. 

In the above-studied two-band model of Sec. 2, the 
relaxation time fi2 equals to f of Eq. (12), and Pi2(k) = 
ak/X. It should be emphasized that both f„ n < and r n are 
inversely proportional to the density of impurities and to 
the matrix elements of scattering from impurities. 

We consider first the terms in Eq. (35) with diago- 
nal elements of the Green functions (i.e., m = w! and 
n = n'). In the static limit of u> — ► 0, the nonvanishing 
contribution comes from the states close to the Fermi sur- 
face, belonging to the same energy branch Ek n (in other 
words, we should also take n = m). Then we find 



I A 



; EE 



g/Cgkn) 

de 



iPx)nn (^y) 



y J nn ' n 



(37) 



This contribution to the Hall conductivity is proportional 
to the relaxation time of electrons. Correspondingly, it 
has the form of the skew scattering mechanismi& 
Let us consider now the contribution from terms with 



n = m and n' 
(36) we find 



in Eq. (35). For u> — > and using 



a IB = e 2 

u xy c 



EE 



df(E kn ) \ , 

T nn * Ilm &n) i^n ~ 4m) 



(38) 



As follows from Eq. (38), this contribution does not de- 
pend on the impurity density and on the scattering am- 
plitude because both inverse relaxation times l/r„ and 
1 /f nn i are proportional to the same impurity density and 
scattering amplitude. Then the ratio T n /f nn i can be pre- 
sented using only the band parameters like in the case 
of model in Sec. 2, see Eq. (33). Hence, the contribution 
(38) should be also identified as an intrinsic mechanism 
even though it was calculated using the impurity scatter- 
ing. In this result we recognize the well-known property 
of the side-jump mechanism^ The only difference from 
the earlier considered models is that the SO interaction 
is included here into the Hamiltonian of free electrons 
instead of the impurity potential. 

Coming back to Eq. (35) with n =/= m and n =/= m , 
we find one more contribution to the AHE, which has 
been discovered recently«A&iZ'i£ and also called intrinsic 
mechanism of the AHE 



e 



« 1_ L m 



k n^m 

x Gk mm (e + w ) (vy) m n G knn (e) , (39) 

where we omit some small corrections related to the off- 
diagonal elements of the Green functions. After inte- 
grating over energy, this result can be presented in the 
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following fori 



,16.17.18 



fdA y (kn) dA x (kn) 



dk„ 



where 



A a (kn) = — i ( kn 



dk a 



(40) 



(41) 



is the gauge potential in the momentum space related to 
the transformation of the Hamiltonian Hq to the diagonal 
form. In a general case, this transformation is local in 
the k-space, leading to the nonvanishing gauge potential 
A a (kn). 

In the model of 2DEG with Rashba Hamiltonian we 
can calculate explicitly the eigenfunctions 



(k,±| = 



'A(fc) ±M 
2\{k) 



1. 



ia(k x — iky) 
M±\(k) 



(42) 



Then using (41) we find the gauge potential 



A(k,±) 



a kn 



a k-r 



2\{k) [M ± X(k)} ' 2A(fc) [M ± X(k)] 



and from Eq. (39) we come again to the same result of 
Eq. (25). Note that (43) can be also found as a gauge po- 
tential corresponding to the local transformation of vec- 
tor field (28) to the homogenous field oriented along axis 
z in the momentum space (like in the case of local trans- 
formations in the real spaced). 

In the 2D case, the flux of curl of the gauge potential 
A(k, ±) in Eq. (40) through the surface -E?k,± = " can be 
presented as the circulation of vector A(k, ±) along the 
circle (Fermi surface in 2D). In other words, the contri- 
bution of the filled states below the Fermi surface can be 
also presented by the integral of the gauge field over the 
Fermi surface. Recently, it was shown by Haldanc 28 that 
such a reduction of the integral in momentum space takes 
place in any dimensionality, in accordance with the Lan- 
dau concept of the Fermi liquid stating that the transport 
properties are fully determined by the properties of elec- 
trons near the Fermi surface. 



in the limit of small impurity density. We demonstrate 
that both contributions can be of the same order of mag- 
nitude. 

As we see from Eqs. (25) and (31), the Hall conductiv- 
ity a xy is proportional to the second order of SO coupling 
a. The even order in a is related to the form of SO inter- 
action in Eq. (1), for which the space inversion, k — > — k 
is equivalent to the change of sign a — > —a. Indeed, the 
nonvanishing integral over momentum in Eq. (6) requires 
the integrand to be invariant with respect to k — > — k, 
which excludes any odd in a factors in o~ xy . 

It is worth to note that the impurity-related mecha- 
nisms of AHE (both side-jump and skew scattering) lead 
to the non- vanishing off-diagonal conductivity in the first 
order of ai The SO interaction included in the impu- 
rity scattering potential has the form of igV^cr (k x k'), 
where g is a constant and Vtk' is the matrix element of 
impurity potential. This type of SO interaction is invari- 
ant with respect to spatial inversion. Thus, in the case 
of small SO interaction, the side-jump or skew scattering 
can be dominating as they both are of the first order of 
SO interaction. 

The essential point of our consideration is that we do 
not neglect the impurity scattering leading to a finite re- 
laxation time of electrons at the Fermi surface. Hence, 
the transition to the " clean limit" in our approach means 
that we should take the limit of small density of impu- 
rities only in the final formula for conductivity after the 
limit of oj — ► 0. The opposite way of calculation would 
lead us to break of continuous transition between the 
" clean limit" and any small impurity density. This point 
was completely neglected in the previous consideration 
of the intrinsic mechanism of AHE. Fortunately, it does 
not affect the "geometric ingredient" of AHE, a™, but, 
as we demonstrated in this paper, it is crucial for the 
competing contribution a xy . 

Note added: After completing the manuscript we were 
informed that similar results have been also obtained 
using a different method by Sinitsyn, Niu, Sinova and 
Nomura^ We thank Jairo Sinova for this information. 
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